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Abstract 

We examine the potentially very promising signal Bg —>■ fJ'^fi^ in supersymmetry 
with large tan/3 in a top-down approach starting from the best fits of an 5'O(10)-like 
model studied recently. Our results go beyond minimal flavour violation investigated 
in previous works. We show that the absolute best fits provide a signal for Bs fi^^^ 
at the borderline of the present limits and hence the ongoing search at the Tevatron 
will start having an impact on the global analysis of this class of SUSY models. We 
discuss the implications of a measurement of Bg —>■ fi^fJ'^ for restricting the parameter 
space of gauginos and sfermion masses, and of signals in other channels B^^s ■ 
We also discuss correlations of Bg /i^/i^ with the CP-odd Higgs mass, sin(/? — a) 
and 6 — s> 57 in S'O(10)-like models. 
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1 Introduction 



Ideas of unification and the origin of flavour have been under investigation for a long 
time and many different models have been proposed in the last twenty years. Yet in 
the diversity of different approaches a class of unification models can be recognised 
which is remarkably simple at the unification scale. We call this class SO(10)-like 
unification models. In these models the effective theory at the unification scale as- 
sumes that the Standard Model (SM) gauge couplings unify to a per cent level, third 
family yukawa couplings are all of order unity and the remaining flavour structure 
originates in a small set of higher-dimensional superpotential operators keeping the 
supersymmetry (SUSY) breaking sector of a model flavour bhnd. We note that actual 
models which fall into this category often assume lower symmetry than SO(IO), e.g. 
models based on the Pati-Salam gauge group or the MSSM gauge group generated 
by a string theory in higher dimensions are often found in this class of models. 

It has been recently pointed out that if the Minimal Supersymmetric Standard 
Model (MSSM) is the effective theory describing nature above the scale 100 GeV and 
tan/3 = / (H^) = Vu/vd is large, a pure leptonic Bg — > decay has a very 

strong case to emerge among the first indirect signals of supersymmetry (SUSY) [1]. 
This is because the decay signal should be very clear at the Tevatron or LHC and 
also because the SM branching fraction is suppressed down to 10"^ while the rate 
can be enhanced when considering SUSY extensions. In particular, this occurs due 
to large couplings of the down-type quarks and charged leptons to the MSSM higgs 
states if tan/3 is large. Thus it is important to analyse this decay in a full SUSY 
theory and not just in terms of the minimal flavour violation that assumes that Vcb is 
the only source of the 32 transition, as has been done in the past. In the full context 
of complete unification models it means that the 32 flavour structure is restricted by 
the fermion mass ratios m^/m,-, ms/rrih and rric/mt, small value of Vcb, large U^s, 
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b ^ branching ratio and possibly other low energy observables and constraints. 
Although these constraints do not determine the 32 sector uniquely they do provide 
for a realistic prediction of observables like the Bg decay rates. 

In this letter we present the results of such a complete top-down investigation 
based on the best fit predictions obtained in a recent global analysis of a complete 
SO(10)-like model [2]. The best fits obtained in this work give a very good agreement 
with the observables related to 32 flavour sector and satisfy all laboratory experi- 
mental constraints on superpartner masses. Here they serve as our starting point 
since they provide us with all the MSSM couphngs at the low-energy scale. Within 
this framework we study the implications of a possible measurement Bg —>■ A*"*"/^". In 
particular, we discuss the related processes Bg — > t'^t~, Bfi — > ii'^iJ,~, B^ t~^t~, 
and show the correlations with Bg — > IJ''^IJ'~ ■ We discuss the implications of a mea- 
surement of Bs A*"*"/^" for restricting the parameter space of gauginos and sfermion 
masses, and also discuss correlations of Bg A*"*"/^" with b ^ and the CP-odd 
Higgs mass. 

After this introduction the letter continues in section 2 with a brief theoretical 
section on the evaluation of Bg — > decay rate in top-down approach. In section 
3 we give some brief discussion of the SO{10)-\ike model we study. Section 4 contains 
our numerical results, and a discussion of the implications of a signal for Bs —>■ A*^A*~ 
mentioned above. Section 5 concludes the paper. 

2 5, -> 

We emphasise that in a top-down approach the tree- level MSSM couplings are de- 
termined from high energy boundary conditions, and do not have to be determined 
by an iterative proceedure as in bottom-up approaches. In particular, in terms of 
effective vertices / and g, which are matrices in flavour space, after heavy sparticles 
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are integrated out the lagrangian can be written down as 



(0) 
R 



+ h.C. 



(1) 



At tree level down-type quarks d^^R only couple to down- type Higgs and / = 
= 0. Yukawa couplings vj^^ can be read out as a straightforward prediction of a 



unified model. and the mass matrix m^^' — Y^'^'vd can then be simultaneously 
diagonalised with eigenvectors d^^\. At one- loop level / and g have to be computed 

and the mass terms relevant for this discussion become ^ 
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using the same basis. Clearly, if v„ S> Vd sizeable corrections to the mass eigenvalues 
[3] and mixing matrices [4] are generated. Furthermore the 3-point functions in (1) 
and mass matrix in (2) cannot be simultaneously diagonalised [5]. If we write eq. (1) 
as 
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then the first bracket of cq. (3) is in a form which is similar to that of the mass 
matrix and therefore is diagonal when c?^'*^ are rotated into corrected mass eigenstates 
d^R — ^d'^^^^d^L^R- This is not true for the last bracket which becomes a source of 
flavour changing, 



-FCNC 



-d 



{1) 
Ri 



Vd 



d^l] + h.c., 



(4) 



It is now explicit that its origin comes from the interaction djiH^*d\, not present at 

tree level. Moreover, the flavour changing couplings get enhanced by an explicit factor 

tan/3 on top of any tan/3 scaling present in g^. In the leading order in tan/3 the g 

matrix can in fact be related in a simple way to the finite non-logarithmic mass matrix 

corrections, gij — (5m^™**'^)ij/vu, computed for the first time in [4]. Due to — 
^ Terms due to wavefunction renormalisation do not contribute to flavour changing. 
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v„+(//V+/i°Ca+iA0c^+iG'%)/V2 and = Vd+{H^Ca-h°Sa+iA^Sf3-iG%)/V2 
we can write 



TTU* TjO 



1 1 

V2 



hPca-p — i A 



(5) 



where = sin a. = cos a, etc. Wc can thus identify effective vertices bnSLH'^, 
hRSihP and hRSL^P involving 6 to s transitions mediated by neutral physical higgs 
states. We note that with large tan /3 the coupling to the pseudoscalar ^4° is always 
large while the CP-even states, and iJ*^, have couplings which depend on the CP- 
even higgs mixing angle a. The Goldstone mode is cancelled in the equation above 
and thus the effective vertex with the Z boson is absent at this level. 

In the MSSM with large tan/? the dominant contribution to Bg comes 
from the penguin diagram where the dilepton pair is produced from a virtual Higgs 
state [1]. After the SUSY partners are integrated out we are left with the effective 
vertices determined above. Thus in combination with the standard tree-level term 
i^UH — —ydR^'LHd'^h.c. the dominant tan (5 enhanced contribution to the branching 
ratio turns out to be 
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V {miooJ 

where matrix 5ml {^L,i?} basis, and is defined by 
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rub is the b quark mass at scale Mz in the effective SU{3)c x U{l)em theory, the 
constants are Mioo = 100 GeV and M^q — 170 GeV and the numerical value is 
obtained from 

tbPbM% 0.0420.03532502 
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(8) 



Modification for other B^. — > decays is trivial. We note that each of these 

branching fractions actually scales down as tan^ /3 for lower values of tan (3: additional 
powers of tan /? enter due to the explicit presence of lepton yukawa coupling and 
mass matrix corrections Sm,'^^"'^^^ / nib (or, equivalently, yukawa coupling y^. in g). 

3 An 5'O(10)-like Model 

Our results are based on the model analysed in [2]. The model was defined be- 
low the SO (10) breaking scale, where the gauge group was broken to its maximal 
Pati-Salam subgroup, and the flavour structure of the model was determined by op- 
erators which respected the Pati-Salam symmetry. Universal gaugino masses M1/2 
and sfermion masses mp were assumed, and we allowed for D-terms and non- universal 
Higgs masses. Throughout this work the trilinear parameter was kept fixed at = 0. 
More details concerning the model can be found in [2] , however the Yukawa matrices 
which enter at the unification scale are listed below for completeness: 
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where A = 0.22 is the Wolfenstein parameter, and a and r are order unity coefficients 
which are precisely determined in the global fit to give excellent agreement with the 
observed quark and lepton masses and mixing angles. The numerical Clebsch factors 
are shown explicitly. Yukawa unification is not exact, with for example dropping 
down to 0.7 for the best fits, although we keep tan/3 = 50 fixed in our analysis. 

The essential features of the fiavour theory clearly include a large off-diagonal 
neutrino Yukawa coupling 1^23 ~ 1, to generate the large atmospheric mixing angle, 
however in the quark sector the Yukawa matrices have small off-diagonal entries, and 
are not required to be symmetric. The ffavour structure of this model is therefore 
typical of many SO{10) models, and has no particularly unusual features, although 
of course wc cannot claim it is generic since each 5'O(10)-like model will differ in the 
details of its flavour structure. 



with no bounds yet established for the r final state processes. Looking to the future, 
the Tevatron will bring us further results for Bg decays with the prospect of a CDF 
bound in the region of Bt{Bs — > < 10~^. By comparison the standard model 
predicts Br{Bs A'.Ai)sM ~ 3.5 x 10~^ [6]. 

To obtain predictions for such processes, we have performed a top-down global 
analysis of the 5'O(10)-like model outlined in the previous section. The analysis 



4 Results 



We first summarise the experimental limits for the processes of interest: 



Br{Bs fifi) < 2.0 X 10-^ [CDF] 
BT{Bd ^ fifi) < 6.1 X 10"^ [Babar] 



(10) 



(9) 
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[2] yields two distinct best fits, which we call Minimum A and Minimum B. The 
higgs spectrum in Minimum B is heavy, mostly above the TeV scale and will not be 
considered in the discussion below. The Higgs spectrum of Minimum A was found to 
be more interesting for our present study with masses at the 100 GeV scale. Hence 
it is the results from the unaltered fits of Minimum A which we present in this letter. 

The numerical results for the processes Bs — > A* and Bg — > t'^t~ are dis- 
played in fig. 1. Similar results for Bd ,t~^t~ are given in fig. 2. These 
results are presented as contour plots in the mp — M1/2 plane with a fixed value of 
fjL — 120 GeV(left panels) and = 300 GeV(right panels). When comparing these 
contours with eq. 6 we find that a significant suppression is obtained from the ratio 
<^"^d 32/^71.6^5. This comes purely from fitting the b quark mass, Vcb and b — > 57. 

The upper two panels of fig. 1 display contours of Br{Bs fifi) with fj, — 120 
and 300 GeV, and show values quite close to the current limits, and well above the 
standard model predictions. The Higgs mediated contribution in the SUSY model 
clearly dominates over the standard model contribution and for fi = 300 GeV, with 
low Ml/2, it can even exceed the present CDF limit. An improved limit of 10~^ would 
be very restricting and could probe Higgs masses into the range, m^o = 150—300 GeV. 
As for the process, B^ jJLji, fig. 2 shows that the present bound is satisfied by both 
II values over the entire displayed plane. 

Inspection of figs. 1 and 2, reveals that the branching ratios for Bg ^ yU^A^", r+r" 
are sensitive to the universal gaugino mass M1/2, but not to the universal sfermion 
mass mp- Inspecting the m^o panels of fig. 3 we see that it has a very similar M1/2, TUp 
dependence. This is exactly as expected with a lighter mediating Higgs leading to 
larger branching ratios. 

The branching ratio for Bg^d 'T't is enhanced by a factor of {yr/ypiY ~ 100 
compared to the muon final state processes, as can be seen in the lower panels of 
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fig. 1 and 2. This makes the tau final state processes very attractive for experimental 
discovery. The difficulty comes with the required detector resolution to measure tau 
decays. If this problem could be solved at future experiments then these tau final 
state processes could become the primary signal for indirect SUSY searches. 

Fig. 3 contains corresponding contours of m^o in the upper panels and the quantity 
sm{P — a), which controls the couphng of the hghtest CP-even Higgs scalar couphng 
to the Z, in the lower panels. The numerical predictions for the best fit point at 
Ml/2 — 450, mp — 500 GeV (indicated by an asterisk in the figures) are given in 
Table 1. 





fi = 120 GeV 


// = 300 GcV 


Ml/2 [GeV] 


450 


450 


niF [GeV] 


500 


500 


Bs 


1.5 X 10-6 


5.9 X 10-6 


Bs^rr 


2.6 X 10"^ 


1 X 10-3 


Bd 


1.5 X 10-^ 


5.8 X 10-^ 


Bd TT 


2.7 X 10-5 


1 X 10-^ 


niAo [GeV] 


102 


102 


sin(/? — a) 


0.22 


0.15 



Table 1: Table of branching ratios for Bg^d IJ-'^ IJ-' , CP-odd pseudoscalar mass 
m^o, and sin(/9 — a) which governs the lightest CP-even scalar coupling to the Z, for 
the best fit point. 

We now turn to the implications of a possible measurement (or an improved 
experimental hmit) of the branching fraction of Bg /jl/jl for 5'O(10)-like models. 
Fig. 4 and 5 show the effect on various quantities of varying the branching ratio for 
Bs —>■ for three fixed points in the mp — M1/2 plane. 

The upper panels of fig. 4 show the variation of Br(i?s jiji) is varied. As 
Br(i?s — > fill) decreases the increases initially slowly and later rapidly. The initial 
slow increase is understood from [2] where it was observed that the value of for 
the best fit points are insensitive to changes of a few GeV in the Higgs spectrum, 
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which imphes an insensitivity to small changes in the branching ratio for Bg — > 
Hence the points which presently exceed the CDF bound can be forced to satisfy it 
with only a small(~ 0.5) increase in x^- But if the bound was to be lowered to 10~^ 
then this would no longer be possible with Ax^ ~ 3. Hence the low M1/2 region of 
the II = 300 GeV plane will be ruled out and the best fit region would move towards 
larger M1/2. 

The lower panels of fig. 4 display the variation of my^o as Bt{Bs nn) is varied. 
As expected m^o increases smoothly as By{Bs /i/i) decreases. Note the strong 
correlation of the CP-odd Higgs mass with Bic{Bs — > /jl/jl), which for a fixed value of 
H is quite insensitive to and M1/2. 

The main contribution to the increase in seen in fig. 4 is due to & — > 57 not 
being fit well. The lower panels of fig. 5 show the variation of Br (6 — > 57) against 
By{Bs — > ii/i) and show a clear correlation. This correlation was to be expected as the 
SUSY contribution to each of these processes involves the 23 mixings in the squark 
mass matrix. These panels also show why b ^ sj is the main contribution to the 
change in as the fit to 6 ^ sj changes from within la to almost 2a. 

The upper panels of fig. 5 show the variation of sin(/3 — a) as By{Bs /i/u) is 
varied. In the low M1/2 region, where By{Bs —>■ /i/J,) is near the current limit, sin(/3— a) 
is small and hence the Z-boson couples predominantly to the heavier CP-even Higgs 
rather than the lighter Higgs However sin(/3 — a) very quickly approaches 
unity as the Bt{Bs —>■ ii/i) decreases, corresponding to the standard model limit where 
the couples like the standard model Higgs boson. 

5 Conclusions 

We have examined the potentially very promising signal Bs A*"*"/^" in supersym- 
metry with large tan /3 ~ 50 in a top-down approach starting from the best fits of an 
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5'O(10)-like model studied recently. Our results go beyond minimal flavour violation 

investigated in previous works. Our results show that the absolute best fits provide 
for the Bs signal at the borderline of the present limits and hence the ongo- 

ing search at the Tevatron will start having an impact on the global analysis of this 
class of SUSY models. 

We have discussed the implications of a measurement (or an improved limit) of 
Bs fi^fi^ for restricting the parameter space of gauginos and sfermion masses, 
and of signals in other channels Bci,s ■ We have also discussed correlations 

of Bg — > ii'^ii" with b ^ and the CP-odd Higgs mass. An improved limit for 
By{Bs —>■ ii/i) of around 10~^ would be very restricting and could probe Higgs masses 
into the range, m^o — 150 — 300 GeV, with the Higgs coupling strength sin(/? — a) 
varying very quickly around this region. The possible non-observation of Bs — > A*"*"//" 
at the levels suggested by our study would by no means rule out S'O(10)-like models. 
In the context of the analysis in [2] this would simply highlight Minimum B, with 
its heavier Higgs spectrum and Br(Ss — > ~ 10~^°, as the favoured solution. 
On the other hand we have seen that an actual observation of Bg fi^fi~ at the 
10~^ level is quite plausibly expected in SUSY S'O(10)-like models, with interesting 
phenomenological and theoretical consequences. 
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Figure 1: Contour plots for the branching ratios of the FCNC processes, Bg M^^^ and 
Bs T~^T~ . Each branching ratio is plotted with two different values of the /i parameter. 
The -k marks the best fit point. 
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Figure 2: Contour plots for the branching ratios of the FCNC processes, — > ijl^ iJi~ and 
Bfi T~^T~ . Each branching ratio is plotted with two different values of the parameter. 
The -k marks the best fit point. 
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Figure 3: The upper two panels contain contours of tlie CP odd Pseudoscalar Higgs mass, 
plotted in the mp — M1/2 plane. The lower panels contain contours of, sin(/3 — a), which 
determines the strength of the Z-boson coupling to, h^, the lighter CP even Higgs. Again 
the plots are displayed at different values of The * marks the best fit point. 
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Figure 4: This figure shows the variation of arid the Pseudoscalar Higgs mass, m^o, as 
the branching ratio of Bg varies from 10~^ down to 10~^. Each of the three curves 

are drawn with fixed values of Mi mp- The vertical dashed line represents the present 
CDF bound on Bg — > 
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Figure 5: This figure shows the variation of sin(/3 — a) and the branching ratio for 
h S7, against Bg /^/i. The vertical dashed line represents the present CDF bound on 
Bg — > iJL^ ^~ . The horizontal lines show the central measured value(solid) of Br(6 sj) 
along with the l(T(dashed) and 2(j(light dashed) regions. 
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